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ABSTRACT 

We present a method that we developed to discern where the optical microvariability (OM) 
in quasars originates: in the accretion disk (related to thermal processes) or in the jet (related 
to non-thermal processes). Analyzing nearly simultaneous observations in three different optical 
bands of continuum emission, we are able to determine the origin of sev eral isolated OM events . 
In particular, our method indicates that from nine events reported by iRarmrez et al. 



three of them are consistent with a thermal origin, three to non-thermal, and three cannot be 
discerned. The implications for the emission models of OM are briefly discussed. 



Subject headings: 
quasars:general 

Introduction 



galaxies: active- galaxies:mndamental parameters 



galaxies:photometry 



Optical microvariability (OM; variations with 
small amplitude in time scales from minutes 
to hours) in quasars can be a powerful tool 
to constrain the energy emission models of ac- 
tive galactic nuclei. Some studies indicate that 
OM does not depend on the radio properties 



of q uasars 
PI 



see 



de Diego et al. 1 11998 



Stali n et al. 



Ramirez et al. 



Gupta k, Joshi 



her eafter 
2005t 



2004 : 

20091 hereafter PII). However 



the mechanism responsible for OM has not been 
characterized. Concerning this topic, variability 
studies in the optical bands acquired a particular 
relevance, because the optical/UV excess is iden- 
tified with the emission from the accretion disk 
' ;.g.. iKrishan fc Wiita I Il994 



Lawrence 20051: 



Pereyra et al. l2006tlLi fc Cao ll2008l : IWilhite et at 
2008). Unfortunately, all the proposed scenarios 
are very complex when the emission from addi- 
tional spectral components is considered. 

In photometric studies, the properties of the 
spectral energy distribution (SED) of quasars 
era characterized by the spectral index, which 
is defined as the slope of a curve in a plane 
log(f v ) versus log{y), and calculated as magni- 
tude differences at different bands fe.g. jMassaro et al 



Il998l iTrevese & Vagnetti I l2002h . During vari- 
ability events, changes in the shape of this 
SED can give us important information about 
the emission processes. For instance, evidence 
from variability studies strongly indicates that 
OM has non-thermal nature in both BL Lac 
objects a nd flat spectrum rad i o quasars (FS 



RQs: e.g.. iD'amicis et al. [2002yVagnctti et al. 



2003: Villata et al. 2004b: Villata et al. 2004a: 



Hu et al. I 12006a , IHu et al. I I2006U IHu et al 



2007), but in the case of FSRQs the contribution 



spectrum of these objects (e.g., 


Gu et al. 20061 


Hu et al. 2006b: Ramirez et al. 


2004). Concern- 



ing quasars, short-term variability (variations with 
amplitudes of tenths of a magnitude and time 
scales of w eeks to months) might h ave thermal 



origin (e.g.. lTrevese fc Vagnetti Il2002h . Neverthe- 



less, all these studies usually employ two bands, 
or only the average of color values is used, losing 
thus information about texture of shortest varia- 



(see Giveon et al. 19991 


Webb & Malkan 


Trevese & Vagnetti 2002 


: Vagnetti et al. 



20031) . 

In order to discern where the OM in quasars 
originates, we assume that emission from the ac- 
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cretion disk must be related to thermal processes, 
while the emission from the jet must be related 
to non-thermal processes. During OM events, the 
spectral component responsible for the variation 
may display characteristic color changes. Then, 
we developed a method with quasi-simultaneous 
observations at three optical bands to analyze 
these color changes that accompanied microvari- 
ability events. This paper is organized as follows: 
in Section [5J we shortly comment on the data 
treatment; in Section |3l we define a spectral vari- 
ability index and how to determine the OM origin; 
in Section 01 we analyze the data; finally, we dis- 
cuss the results in Section [SJ while a summary and 
conclusions are given in Section [5] 

2. Samples selection, observations, and 
data reduction 

Details of the selection criteria, observation 
strategy, and da ta reduction have been exposed 

" (|201dh . These data 



in PI, PII, and Ide D iego 



consist of a sample of 22 core-dominated radio- 
loud quasars (CRLQ) and 22 radio quiet quasars 
(RQQs) observed in diverse epochs. Each ob- 
ject in the RQQ sample paired a CRLQ object 
in brightness and redshift minimizing selection ef- 
fects when properties of microvariability between 
both samples were compared in PII. Here, we will 
only discuss the data corresponding to OM events 
reported in PII. This subsample is listed in Table 
[TJ Four telescopes were used, which are located 
in Mexico and Spain. BVR filters of Johnson- 
Cousins series were used. The observational strat- 
egy consists of monitoring a CRLQ-RQQ couple 
during the same night in overlapping sequences; 
five images of each object in the BVR sequence 
are taken (~ 1 minute of exposure time per im- 
age, although it depends of the object brightness, 
the filters and the telescope used). Each object 
was observed at moderated air masses, always at 
least 30° above the horizon. Standard stars were 
observed to obtain the flux level of the first data 
set of each objects in each night. These stars were 
observed each time that a sequence of each pair of 
objects was complete. This implies that standard 
stars were observed each hour, approxima tely. All 
these stars are taken from lLandolt I (|l992l ). In ad- 
dition, these standard stars are used to perform 
correction by atmospheric extinction. 



3. Quasar SED and SED variability. 

In order to discern where the OM originates, 
the analysis is divided into two main steps: first, a 
SED model is fitted to the initial data set, obtain- 
ing estimations of the values of the parameters of 
each spectral component; then, the observed spec- 
tral variations are modeled using these values. 

The simplest description of the continuum 
emission for quasars is given by a thermal and non- 
thermal emission mixture (e.g . , iMalkan fc Sargent 



1982: iMalkan fc Moore 1119861 ). Thus, at any time 



t, the flux would be determined by the expression 



fvt — Jvnt + fv 



Tt, 



(1) 



where the sub-index t refers to the time when 
a particular data set was acquired; f v t refers to 
the total flux; f unt to its non-thermal component; 
and f U Tt to its thermal component (hereafter sub- 
indices T and n indicate thermal and non-thermal 
components, respectively). Usually, thermal and 
non-thermal components are associated with the 
accretion disk and with the relativistic jet, respec- 
tively Hereafter, we will adopt this assumption. 

At any time t, the contribution from each 
component to the total flux can be described as 
a v t = fvnt/fvt for the non-thermal component, 
and b yt = fvTt/ fvt for the thermal one (note that 
b v t — 1 — GW )- 

In a naive description, the non-thermal emis- 
sion from the jet can be described by a sim- 
ple power law, while the thermal emission from 
the disk can be modeled by a simple black- 
body. Although a simple blackbody is an over- 
simplification for the accretion disk, this is the 
simplest model that fits our dat a (something simi- 
lar has been made previously b y IMalkan fc Sargent 
19821 IMalkan fc Moore~lll986l: among others). 

When a given object is monitored in a partic- 
ular night, we can determine an initial time, to, 
which corresponds to the first data set of observa- 
tions. When a variation is detected we take the 
difference of fluxes between a given data set ob- 
tained at the time t, and that for the first data 
set; then, we normalize with respect to the ini- 
tial total flux. We consider that the variation is 
provoked by only one of the spectral components 
in expression ([1]), leading to the fact that in the 
difference of fluxes only the variable component 
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survives, i.e. 



fvit Ivita 



f, 



(2) 



In this expression, f vit (i = n,T) represents the 
variable component. Some advantages (which are 
easy to verify) of measuring the flux changes by 
this expression are: data transformation to stan- 
dard photometric system is unnecessary; error 
propagation is diminished, and a larger accuracy is 
reached; correction by interstellar extinction and 
cosmological assumptions (as the ones required in 
order to give a value for H ) are also avoided. It is 
important to stress that, when we fitted the model 
to the initial data set, it is necessary to correct for 
interstellar extinction to the first BVR data set; 
while for the spectral variability description, this 
correction is not required. 

On the other hand, we have to consider that 
the observed flux is related to the emission in the 
sources frame through the expression f v = (1 + 



z)L{v \\ + z])/4:1tD 2 l (see iKembhavi fc Narlikar 
1999). Here /„ represents the observed flux at 



the observed frequency v; L the luminosity of the 
source; z the redshift; and Dl the luminosity dis- 
tance. 

In addition, we can represent the brightness 
amplitude of each component by means of a pair of 
constants: A t for the non-thermal component and 
Pt for the thermal one. The injection of new parti- 
cles with similar energy, for the non-thermal case, 
or changes of size of the region of emission, for the 
thermal case, can provoke amplitude changes. 

We initially fit the observed flux at the first 
data set, obtaining the initial values for the tem- 
perature of the disk, the spectral index of the non- 
thermal component, and the contribution of each 
component: T to , at a , and a„t and 6„t , respec- 
tively. Taking the initial data, we can express for 
the non-thermal component the ratio of fluxes be- 
tween any t and to, i.e., 

£«± = _L I/ a« I/ -a t0 ( 1 + jS ja-at Oj (3) 
Junto **to 

from this expression we obtain 



A t 



(4) 



where m vto refers to the magnitude in the stan- 
dard system, corrected by interstellar extinction 
(corrections to these initial magnitudes were car- 
ried out using values reported by the NED); f v 
is a constant of reference flux, corresponding to 
a star of zero magnitude; and v = vb, vv, Vr, 
i.e., the effective frequency in the B, V, and R 
bands. Something similar can be done for the ther- 
mal component, leading to 



fv 



1 



■Tt 



Pt e W (1+x)1 -l 



(5) 



Concerning the contribution from each compo- 
nent, we can give only ay to or bvt , since we can 
write a v t a and b v t (where v — i>b, vr) in terms of 



them, i.e., 



avto !^ 10 -0A( mvto - mvto )fvnt^^ (g) 



fvnto 



while 



Kt a = by t ^ l(T U ^ m -o - m -o > . (7) 

Jv JVTto 

Thus, evaluating expressions (HJ) and ([5]) at 
t = to and replacing them in Equation ([T]), we can 
model the emission at each frequency for the first 
data set to get values for T to , at , and avt ( see 
Figure [TJ. We will test these values when we will 
describe the spectral variations, generating curves 
on a plane defined by zuv,r vs. wb (see Figure 
2). Comparing data and models on this plane, we 
obtain a spectral variability diagnosis plane. 

3.1. Emission models for quasars on diag- 
nosis plane for spectral variability 

Using expression i.e., a power law for the 
non-thermal component, expression ([2]) can be 
written as 

= a vto (n nt {v + zv) a t- at ° - 1), (8) 

while using expression ((5|), i.e., a blackbody as the 
thermal component, expression ([2]) can be written 
as 



^uTt 



v e {hv v /kTt )(l+z) _ ^ 



-b ut0 (n Tt [— } e(h „/fcT t )(i+*) _ 1 



(9) 
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PKS 1510-089 

non thermal emission 

thermal emission 

total with T = 26 000 K 



total 

T =30 000 K 




Fig. 1. — Combination of two components can fit the 
first data set for PKS 1510-089 to obtain Tt , ctt , and 
o-Vto ■ Dotted line corresponds to the non-thermal compo- 
nent, while dashed line corresponds to the thermal compo- 
nent. Spectral index of the power law for the non-thermal 
component has a value of —1, while the temperature of the 
disk would be in 26, 000 K. Continuum thick line represents 
the mixture of these components. However, these data can 
be fitted by different mixtures. Continuum thin lines rep- 
resent the combination of a power law, with index — 1 and 
black bodies of 20, 000 K and 30, 000 K. 



PKS 1510-089 

- thermal variation 

thermal variation with T = 20 000 K 

- thermal variation with T = 30 000 K 





Fig. 2. — Each curve corresponds to a similar thermal 
change for the fits of Figure [I] The curves correspond to 
different values of T tQ . 



where n T t = Pt/Pt and n nt = A t /A to . In this 
way, changes in amplitude are easily obtained by 
varying the values of these ratios from unity, at to j 
to their final values at t. 

Once we get T to , at , and a v t , changes of T t 
and/or n^t, or at and/or n ntl in expression ((2J, 
generate the trajectories in the vjy,R. '■ plane 
(see Figures [2][H]) ■ Simultaneous changes in T t and 
riT t , or at and n nt can generate diagnosis trajec- 
tories with a returning behavior (see Figure [2]). 
A similar case can happen in an expansion of the 
thermal component with a decrease in its temper- 
ature. 

Unfortunately, for the data presented in PII 
it is not possible to establish a unique set of ini- 
tial parameters, due to the lack of data in more 
bands. Figure Q] shows this problem. In this 
figure, we present three reasonable fits for disk 
temperatures of 26,000 K, 20,000 K and 30,000 
K for the first data set of PKS 1510-089. For 
this reason, we must select the most appropri- 
ate values for a to , T to , and a v t , based on the 
typical values found in the literature for simi- 
lar objects. Some investigations suggest that the 



typical maximum disk temperature, close to the 
central object, should be between 20, 000 K and 
30, 000 K, although temperat ures outside of this 
range are not discarded (e.g.. IMalkan fc Sargent 



1982; IMalkan fc Moore I 11989 : ICzernv fc Elvis 
1987t IPerevra et al. 1 12009) . The spectral index 
of the no n-thermal component mig h t be around 
1 ( e.g.. IMalkan fc Sargent] Il982t lElvis et al. 



1989: iBrown et al. Ill989at iBrown et al. Ill989bh 



On the other hand, for blazars, the non-thermal 
component contributes w ith most of the flux (e.g., 



Malkan fc Sargent I Il982l ) , and at low brightness 



levels, it is expected t hat both compo nents con- 



2006). 



tribute the same (e.g.. iGu et al. 

Nevertheless, regardless the exact value for 
these parameters at t = to, we can discern the 
OM origin analyzing the spectral variations. The 
procedure consists of taking several Tt -at -avt 
sets, trying to reproduce the spectral variations 
with each of them. First, the simplest component 
(a power law) is fitted to the initial data, obtain- 
ing a to . From this fit, we look for an explanation 
to the spectral variation, i.e., performing changes 
of amplitude, of the slope, or both, in the diag- 
nosis plane. If the explanation is satisfactory, for 
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simplicity we consider that the OM event is con- 
sistent with changes in this unique component. 
Otherwise, an extra (thermal) component must 
be added to fit the initial data. 



B > 0,10 

0,05 
0,00 
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of curves, as shown in Figures 7. For this reason, 
the estimations for the values of Tt , at , and avt 
must be considered more as a suggestion than as 
the actual conditions for both disk and jet. 



^ 0,10- 




Fig. 3. — Curves show simultaneous variations of at 
and n n t in the non-thermal component for the data of 
PKS 1510-089. Similar changes for ot of -0.7, -1.0 and 
—2.0 are simulated. Temperatures for thermal component 
are 38, 000 K, 26, 000 K and 19, 000 K, respectively. Contri- 
butions of non-thermal component are 0.75, 0.6 and 0.32, 
respectively. The choice of dt = — 1 does not seem to be 
a decisive factor in the variability source discernment. 



Fig. 4. — Three characterization curves for simultaneous 
variations of Tt and fit are shown for data of PKS 1510- 
089. The three curves correspond to similar changes in 
the thermal component for ett = —0.7, —1.0, and —2.0. 
Initial temperatures are 38, 000 K, 26, 000 K and 19, 000 
K, respectively. It is feasible to explain the variability of 
this quasar with the assumption of changes in the thermal 



component. The choice of ct t 



-1 does not seem to be a 



decisive factor in the variability source discernment. 



When this second (thermal) component is 
added to fit the data, two families of characteri- 
zation curves are generated: one for thermal vari- 
ations (as those showed in Figure [2]) and one for 
non-thermal variations (as those showed in Fig- 
ure [3]) . Contrasting the data with each family of 
curves, we can determine whether a spectral vari- 
ation is consistent with a change in the thermal 
component (i.e., changes of T t and/or nxt) or in 
the non-thermal one (i.e., changes of at and/or 
n nt ). 

After comparing the data with all the curves, 
we find that we can set the initial spectral in- 
dex of the non-thermal component, ott , at ~ — 1, 
without affecting our conclusions. So, when 
we include the thermal component, we have 
only T ta and ay ta as free parameters. Then, 
we take this fit as representative of each family 



We have explored how these assumptions can 
affect the discernment of the OM origin. As it 
can be appreciated in Figures [2][6l depending on 
their origin, the curves show a different behavior. 
With different assumptions of values of T to , at , 
and avt > t ne trajectories for non-thermal varia- 
tions are very similar. For this reason, we believe 
that the choice of a to = — 1 has no influence on 
the description of the spectral variation. Expres- 
sion ([8]) shows, in fact, that avt , at — &t i an d 
n nt have more influence on the behavior of the 
trajectories than ctt . Additionally, it is easy to 
see from expression (jHJ that considering another 
power law as the second component, instead of a 
thermal one, does not allow us to reproduce the 
observed variability. 

To illustrate this procedure, we have taken the 
case of PKS 1510-089. First, a unique non-thermal 
component has been fitted to the initial data. 
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0,0 0,1 



0.2 



-0,1 0,0 0,1 0,2 



Fig. 5. — Variation of the thermal component. Curves 
represent similar changes of Tt and nj-, but with o,y to = 
0.60, 0.85, 0.30 (bvt = 1 — a Vt Q )■ A variation from the non- 
thermal component cannot generate behaviors like those 
observed. On the other hand, regardless the exact values 
for Tt , ctt , and ay tQ , the observed variation is consistent 
with changes of the thermal component. 



Fig. 6. — Variation of the thermal component. Similar 
changes in T and nj-, but at different redshifts. z = 0.361 
corresponds to the redshift of the quasar. It is interesting 
to notice that with the same conditions (same Tt , ctt and 
a vt ), objects with different redshifts will have a different 
appearance from the same variation. Additionally, a curve 
representing a change of at and n„t is shown. 



Since this fit does not describe accurately the vari- 
ation, it is necessary to add the thermal emission 
from the disk (Figure [1]). We then fixed «t at — 1 
and calculated the corresponding value for T to for 
different choices of ay fo . Then, we analyzed the 
variations of the thermal component (see Figure 
©. We have used T to ~ 26,000 K (although we 
also explored the cases with 20, 000 K and 30, 000 
K, solid thin lines in Figure [lj and a non-thermal 
component contribution of ~ 60%. 

Afterward, at was fixed at different values, and 
the corresponding values for T to and avt were cal- 
culated. Figures [3] and H] show spectral variability 
curves for changes of the non-thermal and ther- 
mal component, respectively. Curves correspond 
to at = —1, —0.7, —2. As it can be observed in 
Figure [3] for a variation of the non-thermal compo- 
nent, the range of values used for the parameters 
has poor influence on the behavior of the spectral 
variations. On the other hand, the data have a 
behavior that is closer to the one expected from 
thermal changes (Figure 0}. 

Following a different approach, avt was fixed 
at different values (apy = 0.60, 0.85, 0.30), and the 



corresponding values for T to and at were calcu- 
lated. Figure [5] shows different curves for similar 
changes in T t and tit, and for clarity in the Fig- 
ure, we have included only a curve for changes 
in a t and n n , since the other curves are similar 
(as those shown in Figure [5]) . Again, the spectral 
variation has a behavior that corresponds more to 
what is expected from thermal variations. 

Finally, we have explored the effect of redshift 
on variations (Figure ■ At the same values of 
T to , at , and avt , the same variation produced 
at different redshifts (z = 0.0361,0.8, 1.5) causes 
different observable magnitude changes; but with 
greater effect on thermal changes than on non- 
thermal ones. 

From this exploration, we find that variations 
of the non-thermal component cannot explain the 
data; while these are consistent with thermal vari- 
ations. In such a case, the temperature of the 
thermal component varies at least 8 000 K for the 
best fit. 

In summary, in the case of PKS 1510-089, a 
change of the thermal component is the most plau- 
sible explanation for the observed variation, re- 
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gardless of the exact values for T to , a to , and ayt a ■ 
A similar procedure has been carried out for each 
object reported as variable in PII. In the next sec- 
tion we present the results for each object. 

3.2. Temporal resolution of the spectral 
variability 

A problem that has to be solved when dealing 
with color studies refers to the time delay be- 
tween optical bands (e.g . . iGonzalcz-Pe rez et al. 

Il2003i 



1996; Villata et al. 2000; Papadakis ct al 



Hu et al. 1 2006a; Hu et al. 2006b; P apadakis et al. 



2007; IStalin et al. 1120091 ). For spectral OM stud- 



ies, it is particularly important to avoid any pos- 
sibility of a time lag, because a shorter lag than 
the image acquisition time between bands could 
provoke a false spectral variation. Thus, we must 
observe an object at different bands in sequences 
that reasonably guarantee an interval duration less 
than the phys i cal lag s. For the blazar 0716+714, 
Villata et al. I ( 2000h found that a complete se- 



quence between the V and I bands should no t ex- 
ceed ~ 10 minutes. IPapadakis et al. (2003) cal- 
culated a 0.4 hrs time lag between B a nd I bands 
On 



in 



mos t cases. On the other hand, iHu et al. 
(2006b) established a time lag of 11 minutes be- 
tween the BATC e and m bands. Nevertheless 
there is little knowledge about the minimum tem- 
porary scales for microvariability. So, it is nec- 
essary a criterion for maximum acquisition time 
that help us to verify that, during a variation, each 
BVR sequence was make at the same level of flux. 



crete correlation function (e. 


Edelson & Krolik 


1988 


; iHufnaeel & Bregman 


19921; Villata et al. 


2000 


: 1 Villata et al. 2004al). 


but for our observa- 



tional strategy, we developed a time-lag criterion 
based on the change rate of each OM event. This 
criterion permits to be relatively sure that the 
data sets were acquired when each object was at 
the same flux level. 

We introduce a simultaneity criterion that re- 
lates the rate of brightness change with the time 
of acquisition. A simultaneity index is defined by 
the expression 

_Am_ A . 

S, t = ^ff^, (10) 
2.5e„ t 

where Am refers to the maximum variation am- 
plitude; At i, s refers to the time lasted between 



the observations to obtain Am; A.t cap is the re- 
quired time to take a complete sequence of BVR 
images; and e„t is the standard error in each group 
of five images. Then, it will be considered that a 
BVR set of observations is simultaneous when a 
magnitude difference larger than 2.5e is avoided, 
i.e., when the measured error is smaller than a 
marginal variation. Thus, a set of images taken 
at the time t is simultaneous when each S v t < 1. 
As this criterion is applied on each sequence of ob- 
servations for an object in each night, we have an 
average S v for the total group of observations. We 
will show this average in the discussion for each 
object when none of the S v t is greater than 1. In 
other case, each S v t will be showed separately. 

Although this criterion requires that the varia- 
tions have constant change rates, in the next anal- 
ysis we will see that it is possible to make this 
assumption, fulfilling this condition. 

4. Results 

4.1. Model fits and initial data 

The results presented in this section show 
that it is feasible to distinguish qualitatively be- 
tween thermal and non-thermal origin of some 
microvariation events. In those cases where both 
components can explain the spectral variation, 
for simplicity we suppose that the unique non- 
thermal component is responsible for the varia- 
tion. This ambiguity may be eliminated by in- 
creasing the monitoring time and observing more 
optical bands. 

Figures 17-11 17-111 show the procedure described 
in the last section for the case of PKS 1510-089 
for all the objects. The figures show the fits for a 
unique spectral component for each object (panel 
a). The fits for the thermal and non-thermal 
mixed components are shown in panel b. Panel 
(c) shows a plane wy,R ■ for the case of the 
fit showed in panel (a). Finally, panel (d) shows a 
t&v,R ■ plane for the fit showed in panel (b). 
In all cases, a x 2 test indicates that the fits are ac- 
curate. In panels (c) and (d), the error bars have 
been obtained by error propagation in the calcu- 
lation of w v . Additionally, time labels are shown, 
corresponding to to initial data. In all cases where 
two components have been considered, a value of 
at = —1 has been set for the initial slope of 
the non-thermal component. Contribution to to- 
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Table 1: Log of observations, and parameters values from fitting the model with one and two components 



Object 




Observation 




Monitoring 


Magnitude 


Magnitude 






Tto 


a vt 


Variation 


Name 


z 


Data 


Telescope 


Time 


(V) 


Error 


c 


5a 


(K) 


(%) 


Origin 


3C 281 


0.599 


2000 Mar 01 


EOCA 


02:40 


17.40 


0.03 


-0.03 


0.01 


24 000 


46 


Thermal 


3C 281 




2000 Mar 04 


EOCA 


03:45 


17.40 


0.02 


0.13 


0.03 


26 900 


16 


Non-thermal 


PKS 1510-089 


0.361 


2000 Mar 20 


JKT 


04:22 


16.86 


0.02 


-0.23 


0.03 


26 000 


60 


Thermal 


PKS 0003+15 


0.45 


2001 Oct 11 


Mxl 


03:08 


15.36 


0.02 


0.08 


0.02 


19 200 


36 


Thermal (VB 


MC3 1750+175 


0.507 


1999 Jun 20 


JKT 


02:59 


16.18 


0.02 


0.64 


0.02 


33 700 


26 


Thermal (VB) 


CSO 21 


1.19 


2001 Dec 19 


Mxl 


02:32 


17.02 


0.03 


-0.67 


0.11 


29 500 


78 


Non-thermal 


1628.5+3808 


1.461 


2001 Jun 13 


Mxl 


02:23 


16.78 


0.02 


-1.75 


0.02 


10 000 


65 


Thermal 


US 3472 


0.532 


2001 Dec 20 


Mx2 


03:16 


16.25 


0.01 


0.48 


0.01 


27 900 


29 


Thermal (?) 


Mrk 830 


0.21 


1999 Agou 17* 


JKT 


02:11 


17.55 


0.03 


-3.06 


0.49 




100 


Non-thermal 


Mrk 830 




1999 Agou 18* 


JKT 


02:01 


17.56 


0.03 


-2.51 


0.35 




100 


Non-thermal 


Mrk 830 




2001 Jun 14 


Mxl 


01:08 


17.66 


0.01 


-1.81 


0.04 




100 


Non-thermal 



Log of observations for the microvariability events reported in PII (Columns 1-5); values for the fit to the first data set; and possible origin of the variation. 
Column 1 indicate the object name; Column 2, the redshift; the date of the observations are given in Column 3; the used telescopes are given in Column 
4; and Column 5 shows the duration for each monitoring. Telescopes are: EOCA, Estacion Observacional de Calar Alto (Spain); JKT, Jacobus Keptein 
Telescope (Spain); Mxl, 2.01m telescope of OAN (Mexico); Mx2, 1.5m telescope of OAN (Mexico). Table shows initial values for at when only a power 
law is used to fits the first data set (Columns 6-9). The columns indicate the observed V magnitude (Column 6); error in this magnitude (Column 7); and 
the index of the power law and its error (Columns 8 and 9, respectively). The values of the two components fit to the first data set (at = — 1 in all the 
cases), and the possible origin of the variation are shown in Columns 9-12: blackbody initial temperature (Column 10); initial contribution to total flux of 
non-thermal component (Column 11; contribution in the others bands are related to the contribution in the V band, see the text); and the variation origin, 
thermal or non-thermal (Column 12). NOTE: * Variations are only marginal. 
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Fig. 7-1. — Fits to the first data sets and spectral variability: (a) Brst set of data is fitted by a power law; (b) first set of 
data is fitted by a power law and a blackbody mixture; continuum line represents the sum of thermal emission of accretion disk 
(dashed line), and non-thermal emission of the jet (dotted line), (c) the color change is fitted from the fit shown in panel (a); (d) 
the color change is fitted from the fit shown in panel (b); dashed lines represent thermal variations, while dotted lines represent 
non-thermal variations. Errors were obtained from error propagation when to„ is calculated from data, ti (i = 0,1,2...; to 
corresponds to the beginning of each monitoring) indicates the temporary sequence for the data. 



event showed constant rate, our simultaneity cri- 
terion shows that the spectral change is not due 
to observational artifices. 

2000 March 4. We measure m v = 17.40 ± 0.02 
at the beginning of the monitoring. A power 
law can fit these data (see Figure |7^2|) . but with 
a rather low value for the spectral index a to = 
0.13 ± 0.03, with respect to that of March 1. It 
is possible to explain the variation without ap- 
pealing a second component. In such a case, a 
change of the spectral index of Aa ~ 0.031 is re- 
quired (Figure r7-2c|) . However, according to the 
results of March 1, and the fact that at > 0, 
the emission of this object should contain a ther- 
mal component to explain the color variation for 
that night. Since a similar brightness was detected 
in both nights, we included this second compo- 
nent. The thermal component has T to ~ 26, 900 
K, while a vto = 46% and a to = -1 (Figure [T2b}. 
Fitting the variations with two components, the 
solution is degenerated in the sense that both fac- 
tors, thermal and non-thermal, can account for the 
observed variability (Figure [Z2d|; although the 
best fit corresponds still to the non-thermal com- 
ponent. Thus, an index decrease of Aa ~ 0.02 in 
this component, while its amplitude increased by 
n nt2 ~ 1.08, can generate the observed OM. Such 
non-thermal change implies a rate of variation for 
a of ~ 5.3 x 10~ 3 hr -1 . Therefore, it is possible 



to hint a non detectable variation for the second 
data set, with Aa ~ 0.013 (see Figure l7-2dj) . For a 
thermal variation case, it is necessary a decrease in 
T of - 1 800 K, while n Tt2 ~ 1.09. The simultane- 
ity criterion indicates that during the acquisition 
of each set of images the brightness was the same 
(S B = 0.14, S v = 0.48, and S R = 0.39). 

PKS 1510-089. The r eported magnitu de for 



this quasar is my = 16.3 (|Xie et al. 1 120011 ). On 



2000 March 20, we measured m v = 16.86 ± 0.02 
at the first data set. A power law with at = 
-0.23 ± 0.03 fits these data (Figures rT3a|) . How- 
ever, the spectral variation cannot be explained 
appropriately only with a single non-thermal com- 
ponent (Figure !7-3c[) . Thus, a thermal compo- 
nent with T ta ~ 26,000 K and b vto ~ 40% is 
added (Figure [7-3bj) . With this configuration only 
thermal variations can reproduce the observed be- 
havior (Figure [7-3d|) . Nevertheless, the amplitude 
of the thermal component should increase by a 
rather high factor (riT t > 2.5), while the temper- 
ature should fall 8 000 K. In such a case, the rate 
of variation has been AT ~ 1 800 K hr~ x and 
AriT ~ 0.57/i -1 . The return observed in Figure 
!7-3cl can be caused by a variation which is domi- 
nated first by changes in ut and then by changes 
in temperature. This behavior might have been 
evident if observations had been prolonged a few 
hours. With a constant rate of change (as it seems 
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Fig. 7-2. — Fits and variability: continuation. 
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Fig. 7-3. — Fits and variability: continuation. 



to have happened, see the V-band light curve in 
Figure 2(c) of PII), the simultaneity criterion indi- 
cates that observations were carried out satisfying 
the simultaneity criterion (Sb = 0.09, Sy = 0.22, 
and S R = 0.86). 

PKS 0003+15. At the beginning of the obser- 
vations we measured my — 15.3 6 ± 0.02, brighter 
than my = 15.96 reported by iBarvainis et al. 
(|2005l ). A single non-thermal component fits the 
first data set with a to = 0.08 ±0.02 (Figure Oa]). 
However, this fit cannot describes satisfactory the 
color variation (Figure T7-4cp . and a second compo- 
nent must be added. The new fit has a vto ~ 36% 
and T to ~ 19,200 K (Figures Ob]). Both ther- 
mal and non-thermal changes may account for the 
showed variations. The thermal variation is rep- 
resented by the thick line in Figure !7-4dl With 



a thermal variation, as the BV data suggest (su- 
perior panel of Figure !7-4dj) . temperature could 
have increased in ~ 1 200 K, and amplitude de- 
creased by 15% (riTt 3 ~ 0.85). This implies a rate 
of change of 380 K /ir" 1 in T and 0.05ft," 1 for n T . 
The simultaneity criterion indicates that this spec- 
tral behavior is reliable (Sg = 0.02, Sy = 0.19, 
and Sr = 0.01). However, this change would pro- 
duce a detectable variation in the R band, which 
is not detected. Perhaps these data have not the 
required quality, or our simplified model is insuffi- 
cient for explaining this behavior. 

4.3. Radio quiet quasars 

MC3 1750+175. For this object, we measure 
my = 16.18 ± 0.02 at the beginning of the ob- 
servations, i.e., brighter than the value reported 
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Fig. 7-4. — Fits and variability: continuation. 




bv IVeron-Cettv fc Veronl (|200lh (m v = 16.60). 
A power law with ctt — 0.64 ± 0.02 fits this first 
data set (Figure l7-5a|) . However, a second compo- 
nent is necessary to reproduce the observed spec- 
tral variation (Figure !7-5cl) . A mix of thermal 
and non-thermal components fits the data with 
T to ~ 33,700 K and a vto ~ 27% (Figure EHE} . 
Non-thermal changes cannot explain this varia- 
tion. Nevertheless, changes in thermal component 
only can explain the BV data set (Figure !7-5dl) . 
As it happened in the case of PKS 0003+15, a vari- 
ation should be observed in the R band. A varia- 
tion in the thermal component requires a temper- 
ature fall of ~ 5 600 K, while the amplitude should 
increase by 44% (nj- ~ 1.44) (as it is suggested 
from the data in the superior panel of Figure [7-5d[) . 
In such a case, the simultaneity criterion indicates 
the reliability of the spectral evolution (Sb — 0.01, 



Sv = 0.21, and Sr = 0.10); however, this varia- 
tion is something much more complicated to de- 
scribe with our toy model, as it is easy to see in 
the Figure rr5dl 

CSO 21. For this object, we measure a bright- 
ness of my = 17.02 ± 0.03 at the beginning of 
the observ ations, which is similar t o the value re- 
ported by IVeron-Cettv fc Veron I (|200ll ). mv = 
17.0. This first set of data is fitted with ayt — 
-0.67 ± 0.11 (Figure PT6a| . The variation de- 
tected can be explained by an increase of 0.002 
of the spectral index and a decrease in amplitude 
of - 4% (Figure l7-6c|) . It is possible that the OM 
detected for the V and R bands is also present 
in the B band, but the variation in this band 
was not detected due to observational errors. The 
simultaneity criterion indicates that the level of 
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brightness had not changed during observations 
of a BVR sequence (S B = 0.07, Sy = 0.35, and 
Sr = 0.11). For simplicity, we first suppose that 
this OM event results from a variation in a non- 
thermal unique component. Nevertheless, the pos- 
sibility of a second component is revised. In such 
a case, the temperature of the thermal component 
wo uld be T to ~ 30, 000 K, and a V t ~ 78% (Fig- 
ure !7-6b|) . Including this thermal component the 
variations may be generated by any component 
(Figure T7-6dp . In the case of a thermal change, 
the temperature should fall 400 K, while ampli- 
tude would increase by 10% (riT t ~ 1.10). For 
the case of variations of the non-thermal compo- 
nent, the spectral index stays constant, and the 



amplitude increases to n r , 



1.02. 



1628.5+3808. The magnitude of this quasar 
was my — 16.78 ±0.02 at the beginning of the ob- 
servations, i.e., si milar to the reported br i ghtnes s 



my = 16.80 by IVeron-Cettv h Veron I (|200lD . 
A single non-thermal component fits these data, 
with a to = -1.75 ±0.02 (Figure Eg) . However, 
the spectral variation cannot be reproduced by 
changes in this component (Figures !7-7c[) . Then, 
a thermal component is added to model the data 
wit h T to ~ 10,000 K and a vt[) ~ 65% (Fig- 
ure !7-7b[) . A variation of non-thermal origin can 
explain marginally the spectral evolution. On the 
contrary, a thermal variation can reproduce more 
accurately the observed behavior (Figure [7-7d[) . In 
such a case, T would have fallen by ~ 760 K, 
while riT n ~ 1.48. According to the simultane- 
ity criterion, the spectral variation should be real 



(S B = 0.16, S v = 0.12, and S R = 0.14). 

US 3472. When this object was observed, its 
brightness was my = 1 6.25 ± 0.01, i.e., b r ighter 
tha n that reported by La Franca et al. (1992) 
and IVeron-Cettv fc Veron I (|200lh (mv = 16.99). 



Initial data were fitted by a power law with at = 
0.48 ± 0.01 (Figure OH)- Nevertheless, with this 
arrangement we cannot give an explanation for 
the color variation shown in the Figure [7-8cl The 
mixture of a blackbody with T t() ~ 27, 900 K and 
a non-thermal component with ay to ~ 29% can 
also fit the data (Figure !7-8b[) . While a change 
of the non-thermal component cannot explain the 
observed variation, a thermal variation can do it 
(Figure l7-8djl . However, the changes in T and 
riT should be larger in V than in R. In V, the 
temperature would have changed ~ 2 000 K with 
n T t2 ~ 0.87, while in R the change would be of 
only — 1 200 K and n Tt2 ~ 0.92. According to the 
changes of brightness in the V band, the rate of 
the change would be of ~ 600 K /ir _1 . The ob- 
servations between R and B bands were obtained 
with a difference of ~ 8 minutes; then, the tem- 
perature difference would be ~ 80 K. Considering 
the first and the second data sets: Sm — 1.12 
and Sfit 1 = 1.26. In such a case, the temperature 
discrepancy would be explained. The simultane- 
ity criterion indicates that observations were not 
performed at the same level of brightness (con- 
sidering the first and last data sets: Sb = 0.05, 
Sy = 0.53 and Sr — 0.96), because the variation 
in R is very quick, leading to a spurious spectral 
variation. In spite of this lack of simultaneity, the 
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Fig. 7-7. — Fits and variability: continuation. 




spectral microvariation possesses thermal charac- 
teristics. This is evident taking into account that 
the B and V bands data complies with the simul- 
taneity criterion. 

Mrk 830. This quasar showed variations in 
three different nights. Its brightness was similar 
during the three nights (my ~ 17.6), and near 
to the brightness reported by IStepanian et al 



(2001) (m v = 17.29±0.03) and bv lChavushvan et al 



1995) (my = 17.62). A unique component 



marginally fits to the data of this object (Fig- 
ures FT9al I7-I0al 17-llaj) . Although a thermal 
component can be added with T ~ 10, 000 K, 
contrary to the case of 1628.5+3808, the non- 
thermal component is sufficient to explain the 
spectral variation. Then, the addition of a second 
component is avoided. 



1999 August 17. At the beginning of the night 
the brightness of this quasar was my — 17.55 ± 
0.01. But as the first R data set was lost (see 
PII), then the second group was used to fit the 
model. A single non-thermal component fairly fits 
the data, with a to = -3.06 ± 0.49 (Figure FT9a| . 
This value was used for the initial data of B and 
V bands to predict the R value. It is possible to 
explain a microvariation event as observed with 
changes in this non-thermal component (Figure 
!7-9c|) . The spectral index would have increased in 
9.8 x 10 -4 , between the first and the third data set, 
while n n would have stayed constant. The simul- 
taneity criterion indicates that observations were 
taken at the same level of brightness (Sb = 0.09, 
S v = 0.35, and S R = 0.02). 

1999 August 18. The brightness was at the 
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Mrk 830 (1999/Aug/17) 
a„ = -3.06 ± 0.49 



Mrk 830 (1999/Aug/18) 
a. =-2.51 ±0.35 
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Fig. 7-9. — Fits and variability: continuation. Since the 
variation detected on this object can be explained enough 
using a power law, we avoid to make a second fit. In this 
night the first set of data in the R band was lost. 



same level as for August 17, my = 17.56 ± 0.03. 
The first set of observations at the beginning of 
the session can be fitted by a non-thermal compo- 
nent with a to = -2.51 ± 0.35 (Figure OUa)) . The 
variation can be described by changes in this com- 
ponent. In such a case, the spectral index would 
have decreased by a factor of x 10~ 3 and the am- 
plitude would have increased by a factor of 1% 
(n„ ti - 1.01) (Figure FTTOc)) . The simultaneity 
criterion indicates that the spectral variation is 
reliable (S B = 0.27, S v = 0.39, and S R = 0.15). 
However, although the model predicts that the 
variation could be observed in the B band, it also 
predicts that a variations should be observed also 
in the R band, which is not the case. 

2001 June 14- The brightness during this night, 
my = 17.66 ± 0.01, was a tenth of a magnitude 
weaker with respect to the values measured in 
1999 August. Initial data can be described by a 
power law, with a t(t = -1.81±0.04 (Figure lT^TTa]) . 
The microvariability event can be marginally ex- 
plained with a variation in the spectral index 
of ~ 0.001 while the amplitude stays constant 
(Figure I7-Ilc|) . The simultaneity criterion estab- 
lishes that observations have been taken without a 
change in brightness (Sb = 0.17, Sy = 0.50, and 
Sr = 0.37); therefore, the color change is consid- 
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Fig. 7-10. — Fits and variability: continuation 



ered real. 

Although the continuum emission can be fit- 
ted in all cases by a unique non-thermal compo- 
nent, in many of them the broad band spectral 
variation may require to assume the presence of 
a second component. Additionally, considering a 
unique component we have that a to > 0, which 
indicates that more spectral components must be 
acting. 

5. Discussion 

From nine OM detections reported in PII (the 
two possible variations are excluded), our method 
can determine that three events are consistent 
with a thermal variation, and three with a non- 
thermal origin. The other three variations are un- 
determined, i.e., may be associated with either 
thermal or non-thermal processes (the cases of 
PKS 0003+15, MC3 1750+175, and US 3472), but 
with partial characteristics of thermal variations. 

de Diego etal~l (jl998| ) and 



In agreement with 



Ramirez et al. I (|2009), the results presented in 



Section @] show that microvariability does not de- 
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Fig. 7-11. — Fits and variability: continuation. 



pend on radio properties of quasars. Moreover, 
our study shows that OM related to non-thermal 
processes might be detectable not only in radio- 
loud quasars (RLQ) (e.g., 3C 281/2000/Mar/04) 
but also in RQQ (e.g., Mrk 830/2001/Aug/18), 
and that OM related to thermal processes might 
be detected in both kinds of objects, RQQ (e.g., 
1628.5+3808) and RLQ (e.g., PKS 1510-08). 

There are several explanations to account for 
spectral variability originating either in the jet 
or the accretion disk. On the one hand, mi- 
crovariability can be explained by the jet be- 
cause of its physical conditions, and because time 
scales of OM coincide with jet dynamics (e.g., 



Tingav et al~ll200ll : IWiita l[2006h . If all quasars 



can generate a relativ istic jet, as suggested by the 
observational r eports dBlundell fc Rawlings ifeoOll 
Blundell et aL~ll2003t [Ghisellini et al. 112004 ), it is 



not surprising to find the manifestation of non- 
thermal microvariability in RQQ. However, in the 
case of a jet, brightness temperatures exceeding 
the c ritical temperature, T ~ 10 12 K, a re required 



Tingav et al. 



(e.g., 

ences therein; Fuhrmann et al. 



200 it IWiita 1120061 and refer- 
~ 20081) . Although 



this problem may be overcome if large Doppler 



factors are considered ( 
Fuhrmann et al. 1120081 ). 



e.g., iTingav et al. 112001 

The correlation between 
bands represents an additional problem for the jet 
thesis, because correlated variations between op- 
tic /U V and higher frequen c ies should be observed 



(e.g.. iNandra et al. llgOOOt IPeterson et al. I l2000l 



IShemmer et al. I l200ll: iGaskell Il2006h . whic h is not 



the general rule (e.g., Edelson et al. 2000l ) 



On the other hand, variations of diverse time 
scales , including OM, can be produced in accretion 



scales , including um, can pe produced m accretior 
disks I Webb fe Malkan lliooot ITrevese fe Vagnetti 



2001 



2006 



2003 



Trevese fe Vagnetti I 120021: IPerevra et al. 



Wiita 



2006; 



Gu et al. 



see also Mangalam fe Wiita 1993t ~ Fukue 



2006) 



Some of 
associated with disk 
associated with ther- 



these variations may be 
dynamics, or time scales 
mal and sound cross ing time phe nomena (min- 
utes and hours; e.g., Changes in 
the accretion rate can modify the temperature 
of the disk, altering the flux l evel and produc- 
ing short-term variations ( e.g., Webb fc Malkan 1 
20001 IPerevra et al. Il200l I Wiita II2006I ). Pertur- 



bations on the accretion disk can also produce 
variations with ti me scales of weeks, days, and 
Webb fc Malfaml l200d: IWiita 



even hours (e.g. 
20061 ). Recently. 



Chand et al. ( 2010l ) have sug- 



gested that whether the OM is produced by pro- 
cesses related to the jet, then it is expected a 
smaller FWHM of the emission lines in quasars 
with positive detection than in those without OM 
detection. T hese author fo u nd no such behavior. 
Additionally, iGoval et al~l (|2010h have reported 
that the microvariability detected in radio inter- 
mediate quasars has a duty cycle similar to that 
detect ed in RQQs and RLQs by iRarmrez et al. 



(|2009h and IStalin et al. I (|2004h . All these results 



suggest that variations involving the accretion disk 
play an important role. 



ITrevese fc Vagnetti I ([2002J) developed a method 
to analyze the color changes in variability events. 
They determined that short-term variations in 
quasars may be originated by thermal processes 
related to changes in accretion ra te, and hotspots 
in th e accretion disk (see also Vagnetti et al. I 
2003 ). The same con clusion was obtained by 



Webb fc Malkanl (|2000h for variability of quasars 



and Seyfert galaxies. However, although periodic 
or quasi-periodic variations are attributed to cir- 
cling and auto-occulted hotspots in the disk (e.g . , 
Carrasco et al. I Il985t Ide Diego fc Kidger Ill990h . 
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most variations do not display this characteristic. 
Nevertheless, non-periodicity might result of com- 
bination of hotspots and geometric changes, such 
as if the orbi ts in the disk are not Keplerian (e.g., 



Fukue 



20031) . 

On the other hand, the contribution of each 
component to the total flux would play an im- 
portant role to descr ibe the spectral va riations 



(jVagnetti et al. M2003I and iGu et al. 1120061 arrived 



at a similar conclusion) . We need two components 
to explain several spectral variations. In addi- 
tion, when data are fitted exclusively with a non- 
thermal component, we have several cases where 
a to > (3C 281, PKS 003+15, MC3 1750+175, 
US 3472). Since such a behavi or is strange for th e 



1986), 



average SED of quasars (e.g., Elvis et al. 
this result supports the idea that in general a bet- 
ter fit is obtained when a second component is 
added. Cases like those of PKS 1510-089 indi- 
cate that the largest changes of color are due to 
thermal va riations. If a dis k-jet relationship ex- 



ists fe.g.. iFalcke et al. 1119951 ). more complex vari 



ations can be produced by perturbations in the 



disk that are propagated to the jet (e.g., IWiita 
2006 ). Thus, our res ults are in agreement with 



Vagnetti et al. I (|2003l ). who reported that BL Lac 
objects do not show spectral variability as large as 
quasars, for which the thermal component contri- 
bution is larger. 

Our description of the disk emission by a single 
blackbody may be inadequate to explain obser- 
vations; however, more complicated models still 
produce unsatisfactory explanations even for the 
simplest variations. For example, the standard ac- 
cretion disk mod el predicts time lags between op- 



tical bands (e.g., Lawrence 2005; iPerevra et al. 



20061 : IWiita 1120061 and references therein), that are 
not always observed. 

Because for spectral microvariability, variable 
and non-variable components should be consid- 
ered, some spectral components that have not 
been considered by our model can change the 
value of ay ta . These additional components can 
generate spectral variations not contemplated in 
our model; for instance, the contribution from a 
line in only one of the three bands (maybe, this 
can give an explanation to variations as those 
of 3C 281 (2000 March 4), PKS 0003+15, and 
MC3 1750+175). 

In blazars, there are evidences that the spec- 



tral variations are l arger for low flux levels 
(IVagnetti et al. II2003L [Ramirez et a"i~ll2004l : [Gu et al. 
120061 ). This can be understood if at low flux 
level the thermal and non-thermal components 
have a compa rable contribut i on to the opti- 
cal flux (e.g., iRanhrez et al. I 120041) . For BL 



Lacertae, some OM events have been detected 
in which the object turned bluer when it in- 



creases its brightness (e.g., Nesci et al. 


1998; 


Speziali & Natali Ill998; Clements & Carini 


2001; 


Papadakis et al. 1 2003; Villata et al. 


2004b: 


Hu et al. 2006at). These variations should be 



Racine 1970: Gear et al. 1 


19861 iMassaro et al. 


1998: Vagnetti et al. 2003: 


IVillata et al. 1 2004bl 


Hu et al. 2006a). On the other hand, the spec- 



tral changes in FSRQ seem to be larger than in BL 
Lac objects, because the latter have not a remark- 
able contribution from the thermal co mponent 
(|Gu et al. I [20061: iRami'rez et ai~ll2004l) . Then, 



a difference between the two blazars types may 
arise: the thermal component and its influence on 
spectral changes. 

The time-lag criterion that we developed is 
based on the change rate of each OM event for 
our observational strategy. Although this crite- 
rion permits to be relatively sure that the data 
sets were acquired at the same flux level, it re- 
quires that the variations have constant change 
rates. It should be noticed that the light curves 
with a jump, as that of PKS 1510-089 reported in 
PII, can be in fact the result of continuum and 
gradual variations that show such jump because 
the observational strategy. The odd spectral vari- 
ations observed in PKS 0003+15, MC3 1750+175, 
and US 3472 might be caused by the violation of 
the simultaneity criterion. The emission from ad- 
ditional components would also have unexpected 
repercussions on the color changes. 

The importance of orientation on variations 
originated in the jet is not discarded, but processes 
such as those taking place in the disk should be 
considered. Thermal component becomes indis- 
pensable to generate the observed variations, even 
in blazars (particularly at low activity), because 
the accretion disk in these objects must be face- 
on. 

Although variations presented in Section|4]were 
classified as microvariability, we cannot rule out 
that they can be part of a variation of a trend 
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with a larger characteristic time. For example, an 
alteration in the temperature of the accretion disk 
can show a variation with time scales of days, or 
even hours, which might presents a trend that we 
can observe during one night, and interpreted it 
as microvariability. 

6. Summary and conclusions 

The opti cal microvariations in RQQ and RLQ 
reported by iRamirez et al~ ( 20091 ) were analyzed 



tributions on the type of quasar that we observe, 
and if this can be expended to the blazar class. 
This result agrees with the discovery of small 



to determine whether their origin is thermal or 
non-thermal. We suppose that while thermal 
variations are related to the accretion disk, non- 
thermal variations are related to the relativistic 
jet. 

With such a purpose, we developed a method 
to discriminate between color changes by thermal 
and non-thermal processes in quasars during an 
OM event. This method consists in modeling the 
optical broadband continuum of quasars consid- 
ering thermal and non-thermal components, and 
comparing the observed color changes with spec- 
tral variations derived from this model. The main 
result of this work is the possibility to distinguish 
between thermal and non-thermal origin of an op- 
tical microvariation event in quasars. 

We identified that some events are consistent 
with a thermal origin, while some others with a 
non-thermal one. In other cases either component 
can give an explanation to the detected variation. 
Thus, another important contribution of this re- 
search is that, analyzing the spectral microvari- 
ability, we have found that the OM might be gen- 
erated in either the disk or the jet, regardless of 
the radio classification of the quasars, CRLQ or 
RQQ. Additionally, some variations in RLQs are 
consistent with thermal OM and vice versa. Thus, 
our results indicate that microvariability in both 
quasar type might be originated under similar con- 
ditions. 

In addition, although the continuum emission 
can be fitted in all cases by a unique non-thermal 
component, the broad band spectral variation may 
require to assume the presence of a second compo- 
nent, of thermal origin. The relative contribution 
of each component is an important parameter that 
should be taken into account for the description 
of spectral variability. Additional studies must be 
carried on to investigate the effects of these con- 



relativistic jets in RQQs (IBlundell fc R awlings 



2001) an d explains prev i ous re sults (|de Diego et al 



1998 ; iRamfrez et al 



l2009h . Although this 
method yields interesting results applied to our 
data, it is desirable to have a monitoring with a 
better temporal resolution, as well as to observe 
in more bands. More realistic models are also 
necessary. 
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